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The decay of axisymmetric 
homogeneous turbulence 

By J. R. Chasnov 1 


The decay of a homogeneous turbulence generated by an axisymmetric distribution 
of random impulsive forces acting at the initial instant is considered. The impulsive 
forces may be either parallel or perpendicular to the symmetry axis. For impulsive 
forces which result in a k A low wavenumber energy spectrum of the turbulence, it 
is determined that the flow approaches isotropy on all scales of motion at long- 
times provided the Reynolds number is large. However, for the type of impulsive 
forces originally proposed by SafFman [J. Fluid Mech 27 , 581 (1967); Phys. Fluids 
10, 1349 (1967)] in which a k 2 low wavenumber energy spectrum is produced, the 
turbulence approaches isotropy only at the smallest scales and remains significantly 
anisotropic at the largest and energy-containing scales. Nevertheless, a similarity 
state of the flow field establishes itself asymptotically, in which the kinetic energy 
per unit mass of the turbulence decays as f 6 / 5 . 


1. Introduction 

An asymptotic similarity state of decaying isotropic turbulence was predicted 
by Kolmogorov (1941) based on a supposed dynamical invariant of the flow field 
(Loitsianski, 1939). It was later shown by theoretical arguments (Batchelor and 
Proudman, 1956) that the Loitsianski integral is in fact not invariant, and under 
certain conditions of turbulence generation may even diverge (SafFman, 1967a). 
For this SafFman-type flow, a new invariant was discovered and a similarity state 
of decaying homogeneous turbulence at high Reynolds numbers was postulated 
based on this invariant (SafFman, 1967b). Recent large-eddy simulations of isotropic 
turbulence (Chasnov, 1994) have confirmed the existence of this exact similarity 
state to within a few percent. Here, we extend these recent large-eddy simulations 
of isotropic turbulence to a decaying statistically axisymmetric flow field. Our 
main objective here is to determine if a long-time, high Reynolds number similarity 
state of decaying axisymmetric turbulence occurs. Of course, if the axisymmetric 
turbulence approaches isotropy asymptotically, then the earlier found similarity 
state would necessarily be recovered. 

However, a complete return- to-isot ropy of the turbulence is impossible for the 
SafFman-type of flow since the large-scale structure of the turbulence is preserved for 
all times. For instance, if the initial flow field was anisotropic due to some anisotropy 
in the method of turbulence generation, then this initial anisotropy would remain 
in the large-scales of the flow for all times. Of course, the small scales of the flow 
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may become isotropic, but the main scales of interest to us here with regards to 
an asymptotic similarity state are those which contain most of the energy of the 
turbulence; these scales may be more directly affected by the presence of large-scale 
anisotropy. 

Previous closure calculations and numerical simulations have studied the decay of 
an initially axisymmetric turbulence (Herring, 1974; Schumann and Herring, 1976; 
Schumann and Patterson, 1978) in the context of the return-to-isotropy problem. 
The direct numerical simulations performed in the latter two works were necessarily 
limited to low Reynolds numbers, and the computer resources available at those 
times allowed only a resolution of 32 3 . Nevertheless, an approach of the anisotropic 
initial state towards isotropy was indicated. 

This earlier work postulated arbitrary initial states of the axisymmetric turbu- 
lence without considering whether or not the large-scale structure of the flow was 
invariant. Here, following closely the work of SafTman (1967a), we assume that 
our initial flow fields are generated by random axisymmetric impulsive forces at 
the initial instant such as may be generated by laminar flow passing through a 
vigorously shaken grid. We consider impulsive forces which are either parallel or 
perpendicular to the symmetry axis. Both the Saffman-type flow and the original 
Batchelor-Proudman flow will also be considered. In the latter type of flow, the 
large-scale structure of the turbulence is not permanent due to nonlinear transfer 
from small- to- large scales, so that a return-to-isotropy of all scales of motion is 
possible. 

2. Generation by random axisymmetric impulsive forces 

We consider an infinite incompressible fluid initially at rest to which is applied 
random impulsive forces at the initial instant. The continuity and Navier-Stokes 
equations which govern the fluid motion are 

V • u = 0, (1) 

^7 + u • Vu = + i/V 2 u + (S(t), (2) 

Ot Po 

where f is the impulsive force field per unit mass distributed throughout the fluid; 
we assume that this vector force field is a stationary random function of the co- 
ordinate x. Its multiplier <5(t) is the usual Dirac-delta function, which confines 
the impulse to the initial instant. Eq. (2) may be integrated to t = 0 + at which 
time the impulse ceases. Denoting the Fourier transform of u,(x,t) by ib(k,t) and 
the Fourier transform of /i(x) by /j(k), the Fourier components of the resulting 
impulse-generated velocity field at t = 0 + is given by 

Wj(k, 0 + ) = Pij(k)fj(k), (3) 

where -Pjj(k) = 6 t j — kikj] / k 2 , and 6^ is the usual Kronecker-delta. The projection 
operator P»j(k) appearing on the right hand side of (3) occurs because the pressure 
forces respond to the initial impulse so as to maintain the incompressibility of the 
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fluid. The spectral tensor $,j(k, t) of the velocity correlation after the fluid is acted 
on by the impulsive forces is thus given by 


<MM+) = P,m(k)P >n (k)A< mn (k), (4) 

where A4 mn (k) is the spectral tensor of the force correlation, which is assumed 
to exist. The result (4) was previously obtained by Saffman (1967a). We now 
specialize to two particular impulsive force distributions which exhibit statistical 
axi symmetry and could conceivably be generated in a laboratory experiment. 

2.1 Impulse parallel to the symmetry axis 

The physical situation we envision is that of grid-generated turbulence in which 
the grid is randomly shaken along the direction of the mean fluid velocity. The 
impulsive force of the grid on the fluid should ideally be a stationary random func- 
tion of time and of the coordinates in the plane perpendicular to the mean velocity. 
Taylor’s hypothesis can then be invoked to relate the grid turbulence experiment 
to the decaying homogeneous turbulence considered here. 

We shall, without loss of generality, assume symmetry about the X3 axis. The 
spectral tensor of the force correlation corresponding to the above physical experi- 
ment can be written as 


M»j( k) = 6 a 6 j3 M(k), (5) 

and the corresponding spectral tensor of the velocity correlation generated at the 
initial instant thus becomes 


$»>(k,0 + ) = Pj 3 (k)Pj 3 (k)M(k). (6) 

It will be convenient to write the spectral tensor in terms of two unit vectors which 
are perpendicular to k and thus explicitly satisfy the continuity equation (Herring, 
1974). These unit vectors are defined as 


e (1) (k) 


k x n 
|k x n| ’ 


e (2) (k) 


k x e^)(k) 
|k xe< 1 )(k)|’ 


(7) 


where n = (0,0, 1) is the unit vector along the symmetry axis. Using the relation 
Pjj(k) = e^(k)e^(k) + e[ 2 \k)e^(k), it is easy to show that (6) becomes 


<MM + ) = (l - ||) ej 2) (k)ef(k)M(k). (8) 

2.2 Impulse perpendicular to the symmetry axis 

Here, the grid-generated turbulence experiment we envision is the same as above 
except that now the grid is randomly shaken in the direction perpendicular to the 
mean fluid velocity. To maintain the statistical axisymmetry of the force distribution 
and the resulting flow field, we assume that the random impulsive forces acting on 
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the fluid have no preferred direction in the plane perpendicular to the symmetry 
axis. The spectral tensor of the force correlation may now be written as 

M lj (k) = ±(6 ij -6 l3 6 }3 )M(k), (9) 

where the factor of 1 /2 has been inserted for later convenience. The spectral tensor 
of the velocity correlation at the initial instant thus becomes 

*u(M + ) = \ [^(k) - P, 3 (k)P >3 (k)]M(k), (10) 

which in terms of the e-basis can be written as 

*u(k,0 + ) = ± -h ||ej 2) (k)e5 2> (k) M( k). (11) 

2.3 Specification of the spectral tensor of the force correlation 

The initial spectrum of the homogeneous turbulence flow field immediately after 
generation by the random impulsive forces is now completely specified provided the 
form of the spectrum A/(k) of the force correlation is known. For simplicity, we 
assume that M is a function only of the wavenumber magnitude k , The general 
form we choose for M(k) is 

M ( k ) = (£) '*p [- 5 * (£) ] - < 12 > 

where u Q is the root-mean-square velocity of the fluid immediately after generation 
by both types of impulsive force distributions, k p is the wavenumber at which the 
initial spherically integrated energy spectrum of the turbulence is maximum, and 
s = 2 or 4, corresponding either to a Saffman-type flow (Saffman, 1967) or a 
Batchelor-Proudman-type flow (Batchelor & Proudman, 1956), further details of 
which will be presented in the next Section. The normalization constant a s is given 

by 



3. Possible similarity states 

We consider whether a similarity state of the decaying turbulence develops asymp- 
totically in time provided the Reynolds number of the flow is large. Saffman (1967a) 
has shown that when s = 2 in (12), the form of the spectral tensor, (8) or (11), 
persists for all times near k = 0. It is customary to define the associated dynamical 
invariant of the flow field in terms of the leading-order spectral coefficient of the 
spherically-integrated energy spectrum E(k), By spherically-integrating either (8) 
or (11), using (12), the energy spectrum at the initial instant can be shown to have 
the form 
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E{k,Q+) 



(14) 


so that when s = 2, E(k y t) behaves for all times near k = 0 as E(k,t) ~ 2nk 2 Bo^ 
where 


B 0 = 


q2Up 

4nk* 


(15) 


is a dynamical invariant. Normalization constants in (12) have been chosen so that 
(14) agrees with our previous large-eddy simulations of decaying isotropic turbulence 
(Chasnov, 1994). 

A similarity state of the turbulence which is now based on the invariant Bq leads 
directly to the decay of the mean-square kinetic energy per unit mass as (Saffman, 
1967b) 


<u 2 )<xB$t~i. (16) 

The power-law exponent —6/5 was verified to within a few percent in a decaying 
isotropic turbulence (Chasnov, 1994). Saffman (1967a) has further speculated that 
the decay law given by (16) may be valid for a decaying turbulence which does not 
exhibit isotropic symmetry, and one of the purposes of the present work is to test 
this hypothesis. 

When the initial impulsive force distribution is such that s = 4 in (14), the 
spherically integrated three-dimensional energy spectrum of the turbulence follows 
E(k) 27r & 4 i?2 near k = 0, where now B<i = B^ii) is a function of time, as 

shown analytically by Batchelor and Proudman (1957). Lesieur and collaborators 
(Lesieur, 1990), following earlier work of Kolmogorov (1941) and Comte-Bellot and 
Corrsin (1966), have postulated a similarity state of the flow field based on ^(t), 
which was then confirmed within the framework of two-point closure calculations. 
A numerical calculation of i? 2 (f) was also performed by computing an ensemble 
average of a large number (1024) of large-eddy simulations of decaying isotropic 
turbulence (Chasnov, 1993). From the two-point closure calculations and large- 
eddy simulations, the time-dependence of i?2(t) has been shown to be weak relative 
to the overall turbulence decay (the large-eddy simulations of isotropic turbulence 
determined approximately i? 2 (f) ~ t ° -25 ). The decay of the mean-square kinetic 
energy per unit mass may thus be written as 


(u 2 )c (17) 

where the power-law exponent —10/7 is only approximate (a more precise exponent 
of —1.36 results when the time- dependence of B 2 is taken into account). 

Of special interest to the present work is the approach of the initially axisym- 
metric turbulence to isotropy. As a simple measure of the anisotropy of the energy 
containing scales, we define the parameter 7 to be 
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The value of 7 is easily computed for the flow fields generated at the initial instant: 
we determine that 7 ( 0 ) = 8 or 7 ( 0 ) = 2/9 for impulsive forces which are parallel 
or perpendicular to the symmetry axis, respectively. Non-linear transfer processes 
may reasonably be expected to bring the turbulence closer to statistical isotropy 
(7 = 1 ) at later times. 

We have already noted an interesting difference between flows with s — 2 or 4 in 
( 12 ) which may affect the eventual return-to- isotropy of the flow fields. When s = 2 
in ( 12 ), the permanence of the large-scale structure of the flow precludes the entire 
flow from ever becoming isotropic. However, the parameter 7 is primarily a measure 
of the anisotropy of the energy-containing scales so it is yet unclear whether these 
scales will approach isotropy, or whether the anisotropy present in the largest scales 
will affect the asymptotic state of the energy- containing scales. When s = 4 in ( 12 ), 
the large-scale structure of the turbulence is no longer preserved for all times since 
the low wavenumber coefficients of the spectral tensor are no longer invariant. It is 
plausible, although perhaps not entirely obvious, that non-linear interactions can in 
this case produce an isotropic flow field asymptotically in time for sufficiently large 
Reynolds numbers. 

Large-eddy simulations presented in the next Section will shed further light on 
the existence of an asymptotic similarity state in decaying axisymmetric turbulence, 
and the approach (or non-approach) to isotropy of the flow fields. 




4. Large-eddy simulation 

The large-eddy simulations presented here were performed using a pseudospectral 
code for turbulence in a periodic box of length (Rogallo, 1981), and a spectral 
eddy- viscosity subgrid-scale model (Kraichnan, 1976; Chollet and Lesieur, 1981). 
More details about the large-eddy simulation technique can be found in Chasnov 
(1994). Four 128 3 resolution simulations were performed corresponding to an im- 
pulse parallel or perpendicular to the symmetry axis, and s = 2 or 4 in (12) and 
(14). We also choose uo = l and k p = 50. The relatively large value of k p allows 
a similarity state to develop before the integral scales of the flow grow to a size 
comparable to the periodicity length. 

Particular realizations of the impulsive-force-generated velocity fields are con- 
structed as follows. By virtue of the continuity equation, the Fourier components 
of the velocity field may be projected onto the two unit vectors perpendicular to k 
given by (7): 


«i(k) = Mk)e; i} (k) + Mk)eJ 2) (k). (19) 

When the impulsive forces are parallel to the symmetry axis, the spectral tensor at 
t = 0 + is given by ( 8 ) with M(k) given by (12), and we construct random fields 
which satisfy ( 8 ) by setting 
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FIGURE 1. Time-evolution of the power-law exponent of (u 2 ). The solid lines 
correspond to the results of the four large-eddy simulations and the dashed lines 
correspond to the Saffman and Kolmogorov laws discussed in Section 3. The curves 
are labelled by Bq and B2 corresponding to s = 2 or 4 in (12), respectively; and by I 
and II corresponding to impulsive forces parallel or perpendicular to the symmetry 
axis, respectively. 


= 0, ^ 2 (k)=f(i- 

where 9 is a uniformly distributed random number between 0 and 1, chosen inde- 
pendently for each k subject to the complex conjugate symmetry of the Fourier 
components of the velocity field. Similarly, when the impulsive forces are perpen- 
dicular to the symmetry axis, the relevant equations are (11) and (12), and the 
random fields are chosen as 



H 

k 2 


) 


M(k) 


1 i 


exp i2n6, 


( 20 ) 


<M k ) 



exp i2x9i , 


MV = J 



exp 12x02 > 


( 21 ) 


with 8 i and 82 random numbers as above. The three components of the velocity 
field may then be determined directly from (20) and (21) using (19). The particular 
realization of the velocity field constructed above is a generalization of the method 
proposed by Rogallo (1981) for isotropic turbulence. 

In Fig. 1, we present the time-evolution of the power-law exponent (logarithmic 
derivative) of (u 2 ). Time is normalized in terms of the initial large-eddy turnover 
time To, where To = Lo/uo , and Lq is the initial spherically- averaged integral scale 
of the flow, given by Lq = y/i r/k p for s = 2 and Lq = 2y/2n/3k p for 5 = 4 . 
The curves are labeled by Bq or i? 2 , the leading-order spectral coefficients of the 
energy spectrum corresponding to s = 2 or 4 in (12), and the Roman number I or 
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FIGURE 2 . Time-evolution of the anisotropy factor 7, defined in (18). The dashed 
and solid lines correspond to 3 = 2 and 3 ~ 4 in (12), respectively; 7 = 1 for 
an isotropic turbulence. Curves above or below 7 = 1 correspond to impulsive 
forces parallel or perpendicular to the symmetry axis, respectively, (a) Short time 
evolution to t/ro = 100; (b) Long time evolution to t/r 0 = 5000. 

II, corresponding to whether the impulsive forces are parallel or perpendicular to 
the symmetry axis, respectively. The dashed lines are the exact and approximate 
results discussed in Section 3. It is immediately apparent that similarity states in 
good agreement with either the Saffman and Kolmogorov laws develop in all of the 
simulations. 

An evolution of the anisotropy parameter 7(t) as defined in Eq. (18), for short- and 
long-time evolutions, is shown in Fig. 2. Over short-time evolutions, it is apparent 
that all four of the flow fields rapidly become more isotropic than the initial state. 
However, consideration of the long-time asymptotics of 7 clearly demonstrates that 
only when s = 4 (solid-lines) does 7 approach its isotropic value of unity, whereas 
when s = 2 (dashed- lines), 7 deviates substantially from unity at the latest times. 
In fact, the asymptotic values of 7 for these Saffman- type flows are approximately 
7 — ► 1.5 or 7 — ► 0.8 for impulsive forces parallel or perpendicular to the symmetry 
axis, respectively. 

The anisotropy of turbulent eddies should depend on their length scales. To 
observe this effect, we apply the projection of the velocity field defined in (19) to 
construct the following spherically-integrated spectral functions: 


Fj(M) = 4** 2 MM)*i(M)>, Ji(M) = 4 ttP<MM)^(M)), (22 ) 

where * denotes the complex conjugate, and the angular brackets used here denote 
an average over a spherical shell (of unit thickness) in wave space. The usual 
energy spectrum is obtained from E(k,t) = |(i r i(fc,t) + i*2(k,f))- In an isotropic 
turbulence F\(Jfc,t) = i^A:,*), and this will serve as an indication of isotropy for 
scales characterized by a wavenumber magnitude k . In Figs. 3(a) and 3(b) we 
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FIGURE 3. Time-evolution of the spectra defined in (22), for s = 2 in (12) and 
impulsive forces parallel to the symmetry axis, (a) F\(k,t); (b) ^(fc,*). 

plot the F\(k,t) and ^(A;,^) spectra computed from the flow simulation when the 
impulsive forces are parallel to the symmetry axis and s = 2in(12). At the 
initial instant, Fi(fc,0 + ) = 0, as seen from (20). The spectrum of F\(k,t) develops 
rapidly in time and a k 4 low wavenumber spectral form becomes apparent at small 
wavenumbers due to the nonlinear transfer of energy from small-to-large scales. The 
form of the F2(k,t) spectrum as k — * 0 is exactly 4 nBok 2 for all times. 

Using the invariant Bq, asymptotic similarity states can be constructed for the 
spectra F\(k,t) and F2(k,t): 
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FIGURE 4. Time-evolution of the similarity spectra defined in (23), for s = 2 in 
(12) and impulsive forces parallel to the symmetry axis. 


F 1 (it,0 = 5|r«F 1 W, F 2 (k,t) = B$t-*h(k), k = B$ttk. (23) 

The last nine spectra of Figs. 3a and 3b are replotted on the same graph in Fig. 4 
using the scalings given by (23). An excellent collapse of the F\(k) and F 2 (k) 
spectra at different times is observed. Furthermore, the collapse of Fi(fc) and F 2 (k) 
together at large values of k indicates an approach of the small-scale turbulence to 
isotropy. Clearly, however, the flow remains anisotropic at the largest and energy- 
containing scales, as already indicated by the deviation of 7 from unity at large 
times. Apparently, the local value of 7 = 8 for the largest scales of the flow, and 
7 = 1 for the smallest scales results in a balance in the energy containing scales of 
7 « 1.5. 

An interesting result can be obtained from the similarity state given in (23) and 
confirmed by the results of Fig. 4. The form of the F\(k,t) spectrum near k = 0 
can be written as 

F!(k,t)~C(t)k\ (24) 

where, as we have already noted, the k 4 spectrum arises from the nonlinear transfer 
from small- to-large scales. Substitution of this low wavenumber form of the Fi(k, t) 
spectrum into (23), and using the independence of F\(k) on t and Bo , yields the 
exact result 


C(f)oc£ 0 M. 


(25) 
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FIGURE 5. Time-evolution of the spectra defined in (22), for s = 2 in (12) and 
impulsive forces perpendicular to the symmetry axis, (a) Fi(k,t)\ (b) F2(k,t). 

Hence, we have determined the large-scale structure of the turbulence characterized 
by the spectrum F\(k , t), even though this structure is a consequence of complicated 
non-linear interactions. 

Similar ideas apply when the impulsive forces are perpendicular to the symmetry 
axis. In Figs. 5, the time-evolution of the F\(k,t) and ^(A:,^) spectra are plotted. 
The forms of the Fi(k,t) and F2(k y t) spectra near k — 0 are Fi(k y t) ~ 37rJ5oA: 2 and 
*UM) nBok 2 for all times. The rescaled spectra are plotted in Fig. 6 and again 
we find that the largest scales and the energy- containing scales of the flow remain 
anisotropic at the latest times, and the smallest scales approach isotropy. 
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FIGURE 6. Time-evolution of the similarity spectra defined in (23), for s = 2 in 
(12) and impulsive forces perpendicular to the symmetry axis. 


When s = 4 in (12) (and (14)), corresponding to a Batchelor-Proudman- type flow, 
the turbulence becomes isotropic asymptotically and the present results reduce to 
the decay of an isotropic turbulence, whose asymptotic similarity state has already 
been discussed in sufficient detail by Chasnov (1994). 


5* Conclusions 

We have thus demonstrated the existence of asymptotic similarity states of a 
decaying axisymmetric turbulence at high Reynolds numbers. A complete return- 
to-isotropy of the turbulence occurs in Batchelor-Proudman type flows but does not 
occur in Saffman-type flows due to the permanence of the anisotropic structure of 
the large scales. Simulations were performed based on initial flow fields which were 
generated by axisymmetric random impulsive forces acting at the initial instant. 
Such a turbulence generation mechanism has an analogy in active grid turbulence 
experiments, and it would be of interest to see if the results obtained here may be 
reproducible in the laboratory. 


Acknowledgments 

I wish to thank R. Rogallo and A. Wray for allowing me use of their homogeneous 
turbulence code for the Intel iPSC/860 hypercube machine at NAS. I would also like 
to thank C. Cambon, P. Durbin, N. Mansour, and K. Squires for useful discussions 
during the Summer Program. 




The decay of axisymmetric homogeneous turbulence 
REFERENCES 


219 


BATCHELOR, G. K. and PROUDMAN, I. 1956 The large-scale structure of homo- 
geneous turbulence. Philos . Trans . R. Soc. London . 248, 369. 

Chasnov, J. R. 1993 Computation of the Loitsianski integral in decaying isotropic 
turbulence. Phys. Fluids A . 5, 2579. 

CHASNOV, J. R. 1994 Similarity states of passive scalar transport in isotropic 
turbulence. Phys. Fluids. 6, 1036. 

CHOLLET, J. P., AND LESIEUR, M. 1981 Parameterization of small scales of three 
dimensional isotropic turbulence utilizing spectral closures. J . Atmos. Sci. 38, 
2747. 

Comte-Bellot, G. AND CORRSIN, S. 1966 The use of a contraction to improve 
the isotropy of grid-generated turbulence. J. Fluid Mech. 25, 657. 

HERRING, J. R. 1974 Approach of axisymmetric turbulence to isotropy. Phys. 
Fluids. 17, 859. 

KOLMOGOROV, A. N. 1941 On degeneration of isotropic turbulence in an incom- 
pressible viscous liquid. Dokl. Akad. Nauk. SSSR . 31, 538. 

KRAICHNAN, R. H. 1976 Eddy viscosity in two and three dimensions. J. Atmos. 
Sci. 33, 1521. 

LESIEUR, M. 1990 Turbulence in Fluids. Kluwer Academic Publishers. 

LOITSIANSKI, L. G. 1939 Some basic laws for isotropic turbulent flow. Trudy 
Tsentr. Aero.-Giedrodin. Inst. 440, 31. 

ROGALLO, R. S. 1981 Numerical experiments in homogeneous turbulence. NASA 
TM 81315. 

SAFFMAN, P. G. 1967a The large-scale structure of homogeneous turbulence. J. 
Fluid Mech . 27, 581. 

SAFFMAN, P. G. 1967b Note on decay of homogeneous turbulence. Phys. Fluids. 
10, 1349. 

SCHUMANN, U. AND Herring, J. R. 1976 Axisymmetric homogeneous turbu- 
lence: a comparison of direct spectral simulations with the direct-interaction 
approximation. J. Fluid Mech. 76, 755. 

Schumann, U. and Patterson, G. S. 1978 Numerical study of the return of 
axisymmetric turbulence to isotropy. J. Fluid Mech. 88, 711. 




